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The histone deacetylases HDACl and HDAC2 remove acetyl 
moieties from lysine residues of histones and other proteins 
and are important regulators of gene expression. By delet- 
ing different combinations of Hdacl and Hdac2 alleles in the 
epidermis, we reveal a dosage-dependent effect of HDACl/ 
HDAC2 activity on epidermal proliferation and differentia- 
tion. Conditional ablation of either HDACl or HDAC2 in the 
epidermis leads to no obvious phenotype due to compensa- 
tion by the upregulated paralogue. Strikingly, deletion of a 
single Hdac2 allele in HDACl knockout mice results in 
severe epidermal defects, including alopecia, hyperkerato- 
sis, hyperproliferation and spontaneous tumour formation. 
These mice display impaired Sin3A co-repressor complex 
function, increased levels of c-Myc protein, p53 expression 
and apoptosis in hair follicles (HFs) and misregulation of HF 
bulge stem cells. Surprisingly, ablation of HDACl but not 
HDAC2 in a skin tumour model leads to accelerated tumour 
development. Our data reveal a crucial function of HDACl/ 
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Introduction 

The skin is the largest organ of the mammalian body, forming 
a physical barrier against the environment. The external and 
stratified part, the epidermis, is separated from the under- 
lying dermis by the basement membrane and protects the 
body from dehydration and external impacts. The epidermal 
epithelium is estabhshed by differentiation of multipotent 
stem cells (SCs) into the interfollicular epidermis (IFE], the 
hair follicle (HF) and the sebaceous gland (SG) lineages. The 
HFs produce hair in a repeating cycle of hair growth (ana- 
gen) , hair regression (catagen) and a resting phase (telogen) . 
Lubrication of the skin surface is carried out by the SG, 
secreting a lipid-rich sebum. The IFE has a high regenerative 
potential, relying on the constant renewal of its layers — the 
cells are able to self-renew and give rise to the differentiating 
cell progeny. The bulge area, located at the HF beneath the SG 
harbours quiescent SCs, possessing the ability to generate all 
epidermal lineages (Blanpain and Fuchs, 2009). Bulge SCs 
also contribute to the hair cycle and to the wound repair, but 
under normal conditions they do not seem to be involved in 
the homeostasis of the epidermis (Ito et aU 2005). The 
process of terminal differentiation in the IFE is a tightly 
regulated process of gene expression changes ensuring the 
faithful layering of the IFE (Blanpain and Fuchs, 2009] . Cells 
of the basal layer, expressing K5/K14 migrate outward, lose 
their mitotic potential and form the specific keratin markers 
in each layer. Expression of late differentiation markers like 
loricin, filaggrin and involucrin leads to the terminal 
differentiation including apoptosis and formation of the 
cornified envelope. 

As a self-renewing tissue the epidermis is an excellent 
model system for studying the coordinated regulation of 
prohferation and differentiation (Blanpain and Fuchs, 2009]. 
Deregulation of the pathways controlling the transition from 
proliferation to differentiation in the epidermis can lead to 
pathological responses such as skin tumours. Given that non- 
melanoma skin cancers such as basal cell carcinoma and 
squamous cell carcinoma are the most frequent human 
tumours, it is important to understand the basic molecular 
mechanisms and the contribution of epidermal SCs to skin 
cancerogenesis (Boehnke et al, 2012]. Epidermal Hneage 
commitment and morphogenesis are controlled throughout 
embryonic development and adulthood by an orchestrated 
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interplay between different signalling pathways (Blanpain and 
Fuchs, 2006, 2009). In addition, chromatin modifications have 
been shown to play a central role in epidermal 
morphogenesis (Aarenstrup et al, 2008; Ezhkova et aU 
2009; Shaw and Martin, 2009; Ezhkova et al, 2011; 
Botchkarev et al, 2012; Robertson et al, 2012]. Reversible 
histone acetylation is one of the best-characterized chromatin 
modifications and correlates with opening of local chromatin 
structures and transcriptional activation. Two types of 
enzymes, histone acetyltransf erases and histone 
deacetylases (HDACs), control the dynamic acetylation of 
histones and other non-histone substrates. Eighteen 
mammalian HDACs have been identified and classified into 
four groups based on their homologies to yeast enzymes 
(Witt et al, 2009): Rpd3-like class I (HDACl, 2, 3, 8], Hdal- 
like class II (HDAC4, 5, 6, 7, 9 and 10], class IV (HDACll) 
HDACs and the mechanistically unrelated Sir2-like class III 
sirtuins (SIRTl-7]. Since small molecule inhibitors of HDACs 
induce apoptosis, autophagy, cell-cycle arrest and 
differentiation in tumour cells, HDACs are attractive targets 
for anti-cancer treatment. Nevertheless, the precise role of 
individual HDACs in development and tumorigenesis is not 
fully understood. Thus, a fundamental understanding of the 
individual role of HDACs in specific cell types is therefore an 
essential prerequisite for the application of specific HDAC 
inhibitors as therapeutic drugs. 

The highly homologous class I deacetylases HDACl and 
HDAC2 can homo- and hetero-dimerize and are components 
of the Sin3, NuRD, CoREST and NODE co-repressor com- 
plexes (reviewed in Brunmeir et al, 2009). Loss-of-function 
studies in the mouse suggest partial functional redundancy 
for these enzymes in different cell types and tissues 
(Montgomery et al, 2007; Yamaguchi et al, 2010; Chen et al, 
2011 ; Jacob et al, 2011 ; Ma et al, 2012] . Recently, LeBoeuf et al 
(2011) have shown that HDACl and HDAC2 play a central role 
in embryonic ectoderm development. Simultaneous deletion 
of both enzymes mediated by K14-Cre resulted in perinatal 
lethality and impairment of epidermal stratification 
accompanied by increased levels of p53 and impaired gene 
repression by the epidermal key regulator p63. Given a 
potential tumour maintaining role of HDAC 1/2 through 
association with p63 in squamous cell carcinoma, HDACl 
and HDAC2 might also be promising targets in skin cancer 
treatment (Ramsey et al, 2011). By ablating different 
combinations of Hdacl and Hdac2 alleles in the epidermis 
using K5-Cre, we revealed a specific role of HDACl in 
epidermal development and skin cancer. 

Results 

Loss of HDACl or HDAC2 has no obvious consequences 
for epidermal development 

Analysis of HDACl and HDAC2 expression patterns in 
the epidermis revealed overlapping nuclear localization in 
the basal, suprabasal and differentiating cell layers in the 
murine epidermis and HFs (Supplementary Figures SI A and 
B). To examine potential regulatory functions of HDACl and 
HDAC2 in epidermal development, we crossed mice with 
floxed Hdacl or HdacZ alleles (Yamaguchi et al, 2010] to a 
K5-Cre mouse line (Ramirez et al, 2004]. For simplicity, 
Hdacl^^^ K5-Cre mice and Hdac2^^^ KS-Cre mice are referred 
to as Hdacl^^^^^ and Hdac2^^^^^ mice. In agreement with 
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previously pubhshed data (LeBoeuf et al, 2011], Hdacl^^^^^ 
and Hdac2^^^^^ mice showed complete loss of the respective 
deacetylase, and were viable, fertile and displayed normal 
epidermal and HF development (Supplementary Figures SIC, 
D and J]. In accordance with an unaltered epidermal thick- 
ness, immunofluorescence (IF] staining of mouse back skin 
sections revealed comparable expression of the proliferation 
marker Ki67 in Cre-deficient control, Hdacl^^^^^ and Hdac2^^^^^ 
mice (Supplementary Figures SIE and F]. Remarkably, 24% 
of Hdacl^^^^^ mice developed hyperproliferative scars on 
their tails, while no scars were observed on the tails of 
Hdac2^^^^^ mice (Supplementary Figure SIC] . The scar tissue 
showed a significant increase in prohferation as shown by 
quantification of the Ki67 immunohistochemical (IHC] stain- 
ing (Supplementary Figure SIH]. 

HDAC2 protein levels were up-regulated in the epidermis of 
Hdacl^^^^^ mice and HDACl expression was increased in the 
Hdac2^^^^^ epidermis, while mRNA levels were unchanged, 
indicating a post-transcriptional compensatory mechanism 
(Supplementary Figures SlI and J]. To estimate the contribu- 
tion of these two HDACs to the overall HDAC activity, we 
performed HDAC activity assays with epidermal protein 
extracts and tritium-acetate labelled histones as a substrate. 
These experiments revealed no significant changes in total 
cellular HDAC activity in the absence of either HDACl or 
HDAC2 (Supplementary Figure SIK]. 

These data suggest normal epidermal development in the 
absence of either HDACl or HDAC2 with a higher incidence of 
scar formation on HDACl -deficient epidermis. 

Severe developmental abnormalities in the epidermis 
ofHdad^^^^P Hdac^^^^P mice 

To investigate redundant and non-redundant functions 
of HDACl and HDAC2, we generated mice with combined 
ablations of Hdacl and Hdac2 alleles. As reported by LeBoeuf 
et al (2011], mice with simultaneous deletion of Hdacl and 
Hdac2 {Hdacl^^^^^^ Hdac2^^^^'^] were not viable confirming 
an indispensable role of HDACl and HDAC2 for epidermal 
development (Supplementary Figure SIL]. 

HDAC2-deficient mice with a single Hdacl allele {Hdacl^^^^"^ 
Hdac2^^^^^) and compound heterozygous mice {Hdacl^^^^^ 
Hdac2^^^^^) showed normal development and were indistin- 
guishable from Cre-deficient control mice (Figure lA; 
Supplementary Figure SIC]. Interestingly, deletion of one 
Hdac2 allele in Hdacl^^^^^ mice {Hdacl^^^^^ Hdac2^^^^^) 
displayed severe developmental abnormalities (Figures lA 
and C]. Compared to their control littermates, Hdacl^^^^^ 
Hdac2^^^^^ mice were smaller shortly after birth, gained less 
weight from P3 onward and had reduced weight during adult- 
hood (Figures IB and C]. Hdacl^^^^^ Hdac2^^^^^ mice exhib- 
ited progressive alopecia and adult animals had little body hair, 
shorter whiskers and scaly tail regions (Figures IC and E]. 
Remarkably, older Hdacl^^^^^ Hdac2^^"^^^ mice (>4 months] 
but not Hdacl^^^^"^ or Hdacl^^^^"^ Hdac2^/^^'^ mice sponta- 
neously developed papilloma-like lesions on the skin of the tail, 
back, head or extremities (Figure ID]. IHC analysis showed that 
HDAC2 was expressed in all layers of Hdacl^^^^^ 
Hdac2^^^^^ epidermis indicating that the developmental defects 
of Hdacl^^^^"^ Hdac2^^^^'^ mice were not due to restricted 
HDAC2 expression patterns (Supplementary Figure S2A] . 

Deletion of one Hdacl allele in Hdac2^^^^^ epidermis 
or one Hdac2 allele in Hdacl^^^^^ epidermis led to the 
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Figure 1 Severe phenotype of Hdacl^^^^"^ HdacZ^^^^"^ mice. (A) Pictures of adult wild-type control, Hdacl^^^^"^ HdacZ^^^^"^ and Hdacl^^^^"^ 
Hdac2^^^^^ mice. (B) Weight curve of control (^ = 8) and Hdacl^^^^^ Hdac2^^^^^ in = S] mice from PI to P30. (C) Progressive alopecia of 
Maci^/^^P Hdac2^^ + mice with age. (D) Spontaneous skin tumour formation in older Hdacl^^^^"^ Hdac2^^ + mice. (E) Tails of adult control 
and two individual Hdacl^^^^"^ HdacZ^^^^"^ mice (P185, P185, PllO) are shown. (F) Immunoblot analysis of tail epidermal extracts with 
antibodies against HDACl, HDAC2, and p-actin as a loading control. Immunoblot signals quantified by densitometric scanning are shown 
relative to the p-actin signal. n = 4. (G) Total HDAC activity of epidermal extracts from control, Hdacl^^^^^ Hdac2^^^^^ and Hdacl^^^^^ 
Hdac2^^^^^ mice, shown in %. ^ = 4. P = 0.0022 and 0.5172. (F, G) Data are mean±s.d. Source data for this figure is available on the online 
supplementary information page. 



corresponding reduction in mRNA levels (Supplementary 
Figure S2B}. However, immunoblot analysis revealed only 
small changes in the amounts of HDACl and HDAC2 protein 
in the epidermis of Hdacl^^^^'^ Hdac2^^^^'^ and Hdacl^^^^'^ 
HdacZ^^^^^ mice compared to httermate controls (Figure IF], 
indicating that the regulatory cross-talk between HDACl and 
HDAC2 was also functional in these mice. The total cellular 
HDAC activity in epidermal extracts was significantly reduced 
in Hdaci^/^^P Hdac2^/ + ^P but not in Hdaci^/ + ^P Hdac2^/^^P 
mice (Figure IG]. Accordingly, different histone acetylation 
marks including H3K4, H3K27, H3K56 and H4K8 were spe- 
cifically increased in Hdacl^^^^^ Hdac2^^^^^ epidermis 
(Supplementary Figure S2C]. These data demonstrate that a 



single Hdacl allele can compensate for the loss of HDAC2 in 
the epidermis, whereas a single Hdac2 allele fails to balance 
for HDACl deficiency. 

To ask whether this difference was due to different 
expression levels of HDACl and HDAC2 in keratinocytes, 
immunoblot signals for HDACl and HDAC2 were calibrated 
with recombinant proteins and measured. As shown in 
Supplementary Figure S2D, the expression of HDACl and 
HDAC2 protein was very similar as measured with two 
different antibodies for each enzyme. These findings suggest 
that HDACl and HDAC2 have only partially redundant 
functions with a predominant role of HDACl in epidermal 
development. 
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Hdacl^^^^*' Hdac^^^^P mice display disturbed HF 
development 

One obvious characteristic in the Hdacl^^^^"^ HdacZ^^^^"^ 
mice was the reduced pelage and the progressive loss of the 
remaining hair during adulthood (Figures lA and C). 
Therefore, we examined the hair development in more detail 
in the back skin of these mice during HF morphogenesis 
(P5-P10], telogen (PI 8) and anagen (P35] by haematoxylin 
and eosin (H&E] stainings (Figure 2A]. The HF number of 
control and mutant mice was comparable (Supplementary 
Figure S3 A) . However, already during HF morphogenesis (P5 



and PIO] mutant HFs were shorter and disordered compared 
to Hdacf^^ Hdac2^^^ control mice (Figure 2 A and C). At PI 8 
Hdacl^^^^^HdacZ^^"^^^, HFs were significantly longer and 
failed to enter the telogen phase in a synchronized manner. 
Later, the mutant HF became atrophied, lost their morphol- 
ogy and could not enter the anagen phase at P35 as the HFs of 
wild-type littermates (Figures 2A and B]. In line with the 
impaired ability of the HF to fully develop during morpho- 
genesis we detected an increase in p53 expression and 
apoptosis as detected by cleaved caspase-3 staining in mutant 
HFs at PIO (Figure 2D). Thus, p53-induced apoptosis may 
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Figure 2 Progressive hair loss and deregulated HF development in Hdacl^^^^^ Hdac2^^^^^ mice. (A) H&E staining of back skin sections from 
control and Hdacl^^^^^ Hdac2^^^^^ mice at P5, PIO, P18 and P35. Scale bar: 200 |im. (B) Higher magnification of H&E staining of back skin 
sections from control and Hdacl^^^^"^ Hdac^^^^^^ mice at PIO. Scale bar: 200 jim. (C) Quantification of the hair follicle length of back skin hair 
from control and Hdacl^^^^^ Hdac2^^^^^ littermates at different stages. n = 3; P = 0.0353, 0.0002, 0.0336 and 0.0006. (D) p53 and cleaved 
caspase 3-positive (CC3) cells as a fraction of DAPI-stained cells in the HF and in the IFE. n = 4; P = 0.0045, 0.001, 0.0690 and 0.2060. 
(E) Brightfield image of back skin hair from littermate control and Hdacl^^^^"^ Hdac2^^^^^ mice (P35). Scale bar: 100 |im. (F) Back skin HF 
(anagen, PIO) of control and Hdacl^^^^^ Hdac2^^^^^ littermates was stained with an antibody against GATA3. Nuclei were counterstained with 
Mayer's hemalaun. Scale bar: 50|im. (G) Relative mRNA expression of genes important for hair development. ^ = 3-4, P = 0.001, 0.0003, 
0.0001 and 0.0326. (C, D, G) Data represent means ±s.d. 
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considerably contribute to HF degeneration in Hdacl^ 

mice. Remarkably, in comparison to the HFs the 
Hdacl^^^^'^HdacZ^^^^'^ IFE seemed to be more resistant to 
p53-induced apoptosis (Figure 2D; Supplementary Figures 
S3B and C). Hdacl^^^^"^ HdacZ^^^^"^ hair was also thinner 
and the medulla segmentation was partly disrupted 
(Figure 2E). IHC staining for the transcription factor GATA3, 
which is crucial for epidermal lineage determination and 
differentiation of the IRS (Kaufman et aU 2003), revealed 
loss of the protein in Hdacl^^^^^ Hdac2^^^^^ HFs confirming 
an impaired hair morphogenesis programme (Figure 2F]. 
Lhx2, a master regulator of HF SCs (Folgueras et a/, 2013) 
was also down-regulated in HFs of Hdacl^^"^^^ HdacZ^^^^"^ 
mice (Figure 2G). Microarray gene expression analysis of 
control and Hdacl^^^^"^ Hdac2^^^^^ epidermis revealed that 
the disturbed hair development was accompanied by dereg- 
ulation of several genes implicated in hair formation, HF 
development and hair cycle (Supplementary Table SI). Among 
these, several hair shaft keratins {Krt31-Krt86] and keratin- 
associated proteins {Krtapl-Krtap25) as well as the S100a3 
protein, an important structural component of the hair cuticle 
(Kizawa et al, 1998), were found to be strongly diminished in 
their expression in mutant epidermis (Figure 2G). Similarly, 
expression of the important regulators of hair development 
HoxclS and Msx2 was significantly reduced. In summary, hair 
development was severely impaired due to disturbed HF 
morphogenesis and a failure to properly enter the hair cycle 
linked to increased p53 expression and apoptosis resulting in 
degeneration of Hdacl^^^^^ HdacZ^^^^"^ HFs. 

Hyperkeratosis in Hdad^^^^'' Hdac2^^^^P mice 

H&E staining and IF analysis of Kl, K14 and the cornification 
marker involucrin revealed an enlargement of the basal layer, 
spinous layer and stratum corneum in the epidermis of 
Hdaci^/^^P Hdac2^^^^^ epidermis (Figures 3 A, B and D). 
Small proline rich (Sprr) proteins and late cornified envelope 
proteins (LCEs) are important for the process of cornification 
and contribute to barrier function formation (Martin et al, 
2004). Sprr and LCE protein families are encoded by gene 
clusters, which as part of the epidermal differentiation 
complex (EDC) are often coordinately expressed and 
regulated. Several members of the Sprr and LCE protein 
famihes were found up-regulated in Hdacl^^^^^ Hdac2^^^^^ 
mice by qRT-PCR analyses and gene expression profiling of 
epidermal mRNA (Figure 3C; Supplementary Table SI). 
This up-regulation is consistent with the hyperkeratosis 
observed in Hdacl^^^^^ Hdac2^^^^^ epidermis. The observed 
hyperkeratosis and epidermal hyperproliferation was not 
connected to a defect in barrier function, since both transe- 



pidermal water loss (TEWL) test and toluidine blue dye 
exclusion assay gave no significant difference between 
Hdacl^^^^"^ Hdac2^^^^'^ mice and wild-type littermates. 
Also skin barrier genes were unaffected in the mutant epi- 
dermis (Supplementary Figures S4A-C) . These results show a 
deregulation of differentiation-related genes and hyperkera- 
tosis in the epidermis of Hdacl^^^^"^ Hdac2^^^^'^ mice. 

Increased SGs and IFE hyperproliferation in Hdacl^^^^^ 
Hdac2^^^^P mice 

In order to survey the SO lineage, we stained control and 
Hdacl^^^^^ Hdac2^^^^'^ skin whole mounts with Oil Red 0, 
a histochemical dye specific for lipid-containing cells including 
sebocytes. Hdacl^^^^'^ Hdac2^^^^'^ SGs were significantly 
enlarged and filled the region of the HF shaft, which was 
underdeveloped or entirely absent (Figure 3E). The serine 
protease kallikrein related-peptidase 6 (Klk6/BSSP), which is 
predominantly expressed in the SGs (Frye et al, 2003; 
Komatsu et al, 2005), was up-regulated in the epidermis of 
Hdacl^/^^'' Hdac2^/^^^ mice (Figure 3F). 

To determine whether the observed thickening of the 
epidermis in Hdacl^^^^^ Hdac2^^^^^ mice was due to epi- 
dermal hyperproliferation, we analysed BrdU incorporation 
at P5 and Ki67 expression at P30. Both were found to be up- 
regulated in the IFE of Hdacl^^^^^ Hdac2^^"^^^ mice 
(Figures 3H and I; Supplementary Figure S3D). We further 
detected a significant induction of epithelial mitogen [Epgn] 
and adenosine deaminase [Ada], two genes crucial for 
epithelial morphogenesis and proliferation, in the mutant 
epidermis by qRT-PCR (Figure 3G). Keratin 6 (K6), a marker 
for epidermal hyperproliferation is, under normal conditions, 
predominantly found in the companion layer of HFs and 
becomes expressed in the epidermis only upon stress or 
injury (Mazzalupo et al, 2003). Notably, we found up- 
regulation of K6 as well as other stress-associated genes 
including keratin 16 [Krtie], SlOOaS and S100a9 
(Supplementary Figure S3E; Supplementary Table SI). 
However, FACS analysis for markers of immune cells at P5 
indicated no significant increase in immune infiltrates in the 
Hdacl^^^^^ Hdac2^^^^'^ epidermis excluding that an inflam- 
matory condition could be responsible for the observed 
hyperproliferation (Supplementary Figures S4D and E). 
These data show that haploinsufficiency of Hdac2 in the 
absence of HDACl resulted in increased SGs and hyperproli- 
feration of the IFE. 

Changes in lineage determination in Hdad^^^^^ 
Hdac2^^^^^ epidermis 

Since all three epidermal lineages were affected in their 
development in Hdacl^^^^^ Hdac2^^"^^^ mice, we examined 



Figure 3 Hdacl^^^^^ Hdac2^^^^^ mice display impaired epidermal development and hyperkeratosis. (A) IF analysis of back skin sections of 
adult control, Hdacl^^^^^ Hdac2^^^^'^ and Hdacl^^^^^ Hdac2^^^^^ mice (P35) with antibodies specific for K14 and Kl. Nuclei were 
counterstained with DAPI. Scale bar: 50|im. (B) Quantification of epidermal thickness as shown in (A) in littermate controls and Hdacl^^^^^ 
Hdac2^^^^^ mice. n^3. P = 0.0150. (C) qRT-PCR analysis of mRNA expression of Sprr genes in the epidermis from Hdacl^^^^"^ Hdac2^^^^^ 
and control littermates. n^3. P = 0.0190 and 0.0139. (D) IHC analysis of control and Hdacl^^^^^ Hdac2^^^^^ back skin sections for involucrin 
(P35). The nuclei were counterstained with Mayer's hemalaun. Scale bar: 100 |im. (E) Brightfield images of Oil red 0 and haematoxylin-stained 
skin whole mounts of control and Hdacl^^^^^ Hdac2^^^^^ mice (P90). Dashed line marks SG, straight line marks HF unit. Scale bar: 200 |im. SG 
areas for individual HF were quantified and shown on the right panel. ^ = 30. P = 0.0315. (F) Relative mRNA expression of the kallikrein- 
related peptidase 6. n = 3. P = 0.0165. (G) qRT-PCR analyses of Epgn and Ada in the epidermis of control and Hdacl^^^^^ Hdac2^^^^^ mice. 
n = 3. P = 0.0122 and 0.0250. (H) Quantification of IF signals for BrdU in HF (left panel) and IFE (right panel) of control and Hdaci^^^^P 
Hdac2^^^^^ mice (P5) after 2 h BrdU labelling. n = 4. P = 0.0007 and 0.0068. (I) IHC analysis for the proliferation antigen Ki67 in control and 
Maci^/^'P Hdac2^^^^'P tail skin (P30). Scale bar: 100 |im. Quantification of Ki67-positive cells in tail epidermis is shown on the right. n = 3. 
P = 0.0150. (B-I) Mean value with s.d. is shown. 
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the expression of markers for SCs located in the HF bulge 
iCD34, Keratin 15) and hair sheath iSox9, LgrS, Lgr6). To this 
end, we performed qRT-PCR analyses of the epidermis of 
Hdacl^^^^'^ HdacZ^^^^'^ mice and control littermates. 
Strikingly, most analysed SC genes such as CD34, Keratin 
15, LgrS and Sox9 were significantly down-regulated in adult 
Hdacl^^^^^ HdacZ"^^^^^ epidermis (Figure 4A]. The 



R-spondin receptor LgrS is a marker for multipotent SCs able 
to generate ah lineages of the HF (Schuijers and Clevers, 
2012). To investigate its local expression pattern, we 
performed in situ hybridization experiments. In control 
epidermis, LgrS was located in the HF bulb and ORS, 
whereas LgrS was not detectable in HFs of Hdacl^^^^"^ 
Hdac2^^^^^ mice (Supplementary Figure SSA). On the 
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Figure 4 Epidermal hyperproliferation, deregulated epidermal lineage commitment and c-Myc activation in Hdacl^^^^^ Hdac2^^^^^ mice. 
(A) Relative mRNA expression of epidermal stem-cell markers in tail epidermis from control and Hdacl^^^^^ Hdac2^^^^^ littermates. ^ = 3-4. 
P = 0.0054, 0.0006, 0.0022, 0.0081 and 0.0008. (B) IF co-staining of label-retaining (BrdU + ) bulge cells (CD34) in HF (left panels) and IFE 
(right panels) of control and Hdacl^^^^^ HdacZ^^^^"^ mice (P34) after 72 h BrdU pulse (P6) and 25 days chase. Nuclei were counterstained with 
DAPI. Scale bar: 100 |im. Quantification of BrdU is shown on the right of each staining picture. n = 4. P = 0.0241, 0.0030 and 0.0646. (C) Back 
skin sections of control and Hdacl^^^^"^ HdacZ^^^^^ mice (P30) labelled with an antibody against c-Myc. Nuclei were counterstained with 
DAPI. Scale bar: 50|im. (D) Expression of c-Myc in epidermal tail nuclear cell extracts of control, Hdaci^^ + ^P HdacZ ^^^^p and Hdacl^^^^^ 
Hdac2^^^^^ mice was analysed on immunoblots. Lamin B served as a loading control. Immunoblot signals of control and Hdacl^^^^^ 
HdacZ^^^^"^ mice were quantified by densitometric scanning are shown relative to the Lamin B signal on the right of the panel, n = 3. (A, B, D) 
Error bars indicate s.d. Source data for this figure is available on the online supplementary information page. 



contrary, Lgr6 an SC marker that is postnatally contributing 
to SG and IFE growth (Snippert et al, 2010; Schuijers and 
Clevers, 2012) show^ed increased expression in the mutant 
epidermis (Figure 4A). To examine w^hether there was a 
change in SC fate in Hdacl^^^^^ HdacZ^^^^^ epidermis, we 
performed BrdU pulse-chase experiments to analyse slow 
cycling label retaining cells (LRCs) and their progeny. 
Control and Hdacl^^^^"^ HdacZ^^^^"^ mice (P6) were pulsed 
for 72 h with BrdU, followed by a 25-day chase. Proliferating 
cells lose the BrdU label whereas slow cycling LRCs retain the 
BrdU label. To identify SCs in the bulge, co-stainings for BrdU 
and CD34 were performed. As shown in Figure 4B, we 
detected a reduction in BrdU + /CD34+ cells in the HFs of 
Hdacl^^^^"^ HdacZ^^^^^ mice, which is at least in part caused 
by increased cell death (Figure 2D). At the same time, we 



observed an increase in BrdU ^ cells in the IFE and (albeit not 
significant) in SGs of mutant mice (Figure 4B) . These findings 
reveal profound changes in epidermal lineages infhcted by 
enhanced apoptosis in the HFs, increased proliferation in the 
IFE and/or mobilization and redirection of epidermal SCs 
in Hdacl^^^^^ HdacZ^^^^^ mice. The increased proliferation 
in the IFE suggests that the bulge SCs are able to contribute to 
IFE regeneration. This prompted us to perform wound-heal- 
ing assays by introduction of full-thickness excisional 
wounds in the back skin. Notably, mutant mice display a 
delay in wound healing and developed keratoacanthoma-Uke 
lesions with massive keratin deposition around the wound 
margins (Supplementary Figures S5B and C). 

The phenotype of Hdacl^^^^^ HdacZ^^^^^ mice with 
hyperproliferation of the IFE, increased SG formation 
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and impaired hair development is reminiscent of c-Myc 
overexpressing mice, which display epidermal hyperproli- 
feration and differentiation along the SG and IFE lineages at 
the expense of HF development (Arnold and Watt, 2001; 
Waikel et aU 2001). Therefore, we analysed c-Myc 
expression in the epidermis of control and mutant mice. 
Hdacl^^^^^ Hdac2^^^^^ mice displayed increased c-Myc 
protein levels in the epidermis detectable both by IHC and 
by immunoblot analysis, whereas Hdacl^^^^^ HdacZ^^^^^ 
and control epidermis showed no c-Myc induction 
(Figures 4C and D). Accordingly, the c-Myc protein was 
also strongly expressed in spontaneously formed tumours 
of older mutant animals (Figure ID; Supplementary Figure 
S5D]. Collectively, our results obtained with Hdacl^^^^^ 
Hdac2^^^^^ mice are comparable with data pubhshed for 
c-Myc overexpressing mice which display increased prolifera- 
tion and differentiation along the epidermal and sebaceous 
lineages at the expense of HF differentiation (Arnold and 
Watt, 2001; Waikel et al, 2001). 

Deregulated gene expression in Hdad^^^^^ Hdac2^^^^'' 
epidermis 

Given the different epidermal phenotypes in mice with Hdacl 
disruption in the presence and absence of up-regulated 
HDAC2 levels, we compared gene expression profiles of the 
epidermis of Hdaci^/^'P, Hdacl^^^'^ Hdac2^^^^^ mice and 
the corresponding control httermates. A global expression 
analysis showed that in Hdacl^^^^^ mice only 79 genes were 
deregulated compared to control littermates (P<0.05, at least 
two-fold change in expression) (Figure 5A; Supplementary 
Table S2). In contrast, 3749 genes were deregulated in the 
epidermis of Hdacl^^^^^ Hdac2^^^^^ mice suggesting that up- 
regulated HDAC2 can compensate for most of the lost 
regulatory function of HDACl in Hdacl^^^^^ epidermis 
(Figure 5B; Supplementary Table S2). We categorized the 
up-regulated genes by functional gene ontology (GO) analysis 
as shown in Supplementary Table S3. In agreement with the 
observed phenotype, genes in categories such as cell prolif- 
eration and keratinization were enriched in the group of up- 
regulated genes. In a recent pubHcation, Nascimento et al 
(2011) reported a functional hnk between Sin3A and c-Myc in 
the regulation of epidermal development. In the absence of 
Sin3A in the epidermis, c-Myc was up-regulated resulting in a 
phenotype that is reminiscent of Hdacl^^^^^ HdacZ^^"^^^ 
mice. Consistent with this idea we found that more than 
one third of the genes deregulated in Sin3a^^^^P epidermis 
(Nascimento et aU 2011) were also deregulated in the 
Hdacl^^^^"^ HdacZ^^^^"^ epidermis (Figure 5C). Further- 
more, a subset of these genes (166) have previously been 
shown to be deregulated upon overexpression of c-Myc in the 
epidermis and a significant number of genes are commonly 
deregulated in Hdacl^^^^^ HdacZ^^^^^ and c-Myc overex- 
pressing epidermis (Figure 5C) (Frye et al, 2003). These 
data suggest that an impaired Sin3A co-repressor function 
contributes to the phenotype of Hdacl^^^^^ HdacZ^^^^"^ mice. 

Partially impaired co-repressor function in Hdad^^^^^ 
Hdac2^^^^P epidermis 

Therefore, we next investigated by co-immunoprecipitation 
experiments whether HDAC1/HDAC2 co-repressor activity 
was affected in Hdacl^^^^^ HdacZ^^"^^^ keratinocytes. 
Remarkably, deacetylase activities associated with Sin3A, 



the NuRD complex component MTA2 and CoREST were all 
reduced in Hdacl^^^^^ HdacZ^^^^^ keratinocytes (Figure 5D). 
Immunoblot analysis showed reduced protein levels of Sin3A 
and MTA2 in Hdacl^^^^'^ HdacZ^^^^"^ keratinocytes, whereas 
CoREST expression was not affected (Figures 5E and H). The 
reduction in MTA2 and Sin3A protein levels was not caused 
by reduced mRNA expression but by lower protein stability as 
shown in cycloheximide protein stability assays (Supplementary 
Figures S6A and B). In contrast, Hdacl^^^^^ HdacZ^^^^"^ kera- 
tinocytes showed no significant changes in co-repressor asso- 
ciated deacetylase activities (Figures 5F, G and I) . Thus, loss of 
HDACl in the absence of elevated HDAC2 leads to strongly 
reduced deacetylase activity and partial destabilization of 
HDAC1/HDAC2 repressor complexes. 

Sin3A (and associated proteins) has been previously 
shown to cause deacetylation and destabiUzation of the 
c-Myc protein (Nascimento et al, 2011). Of note, the up- 
regulation of c-Myc protein in the Hdacl^^^^^ HdacZ^^^^^ 
epidermis was not accompanied by increased c-Myc 
mRNA levels suggesting a post-transcriptional regulatory 
mechanism (Supplementary Figure S6C). In accordance 
with this hypothesis, HDAC inhibition with the class 
I-specific deacetylase inhibitor MS-275, which has a prefer- 
ence for HDACl (Bertrand, 2010) led to increased c-Myc 
protein levels and protein stabihty without affecting 
c-Myc mRNA expression (Supplementary Figure S6D). 
These data, together with the findings of Nascimento et al 
(2011), indicate a crucial function of the Sin3A co-repressor 
complex and associated HDAC1/(HDAC2) activity in 
preventing inappropriately high levels of c-Myc protein in 
the epidermis. 

Several members of the differentiation-associated members 
of the Sprr and LCE protein families which had been found to 
be negatively regulated by Sin3A are also up-regulated in the 
Hdaci^/^^P HdacZ^^^^"^ epidermis (Nascimento et al, 2011) 
(Supplementary Table SI). To ask how the impaired Sin3A 
co-repressor function affects the regulation of genes com- 
monly deregulated in Sin3A-deficient and Hdacl^^^^^ 
Hdac2^^^^^ epidermis, we performed chromatin immunopre- 
cipitation (ChIP) assays. The association of HDACl, HDAC2 
and Sin3A and the presence of histone H3K9ac and H4ac 
marks was analysed for up-regulated genes representing the 
EDC [SprrZh, SlOOaS, S100a9, LceSh), SG {Klk6), epidermal 
proliferation markers {Epgn, Ada], and the non-regulated 
control VstmZl with chromatin isolated from epidermis of 
control, Hdacl^^^^^ and Hdacl^^^^^ HdacZ^^^^^ mice. In 
wild-type controls, HDACl, HDAC2 and Sin3A were present 
at the promoter regions of SprrZh, SlOOaS, S100a9, LceSh, 
Klk6, Epgn, Ada, HoxclS and Msx2 (Figure 6; Supplementary 
Figures S7A and B). In Hdacl^^^^"^ epidermis, where these 
genes are not deregulated, increased amounts of HDAC2 were 
associated with these promoters in the absence of HDACl. In 
contrast, in Hdacl^^^^^ Hdac2^^ + epidermis the presence of 
Sin3A and HDAC2 at the SprrZh and Epgn promoters was 
reduced concomitant with an increase in local histone acet- 
ylation and gene expression. On the other hand, the down- 
regulated hair development regulator genes [HoxclS, Msx2) 
showed reduced histone acetylation levels upon loss of 
HDACl in Hdacl^^^^"^ HdacZ^^^^"^ epidermis suggesting a 
potential positive role of the deacetylase for the expression 
of these genes (Supplementary Figure S7B). We conclude that 
increased presence of HDAC2 at these target genes can 
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Figure 5 Gene expression and changes in repressor complex function in Hdacl^^^^"^ Hdac2^^^^^ mice. Agilent microarray gene expression 
analysis otHdacl^^^^^ (^ = 4) [A], Hdacl^^^^^ HdacZ^^^^"^ (^ = 3) (B) and the corresponding control mice. RNAwas isolated from adult tail 
epidermis from three animals of each genotype. In all, 79 annotated genes in Hdacl^^ mice (A) and 3749 genes in Hdacl^^^^^ Hdac2^^^^^ 
mice (B) were > 2-fold deregulated. (C) Overlap of deregulated genes of Hdacl^^^^"^ HdacZ^^^^"^ mice v^ith deregulated genes of SinSa^^^^"^ 
mice (Nascimento et al, 2011). Out of this overlap of 550 genes, 166 genes v\Aere deregulated in c-Myc overexpressing epidermis (Frye et al, 
2003). Overlap of deregulated genes in c-Myc overexpressing (OE) epidermis v^ith deregulated genes of Hdacl^^^^^ Hdac2^^^^^ mice. 
Significance of overlaps was determined by a hypergeometric probability test. (D-G) Sin3A, NuRD and CoREST complexes were 
immunoprecipitated from keratinocyte extracts of control, Hdacl^^^^^ Hdac2^^^^^ (D, E) and Hdacl^^^^^ HdacZ^^^^^ (F, G) mice with 
antibodies against Sin3A, MTA2 and CoREST. (D, F) Complex-associated HDAC activity was measured in input and IPs. Controls are set to 
100%. ^ = 4. (E, G) Immunoprecipitated proteins were analysed on immunoblots with HDACl, HDAC2, Sin3A, MTA2 and CoREST antibodies. 
P-actin served as a loading control. (H, I) Immunoblot signals for Sin3A, MTA2 and CoREST in Hdacl^^^^^ Hdac2^^^^^ (H), Hdacl^^^^^ 
Hdac2^^^^^ (I) and the corresponding control keratinocytes were quantified by densitometric scanning and are shown relative to the p-actin 
signal. The values of wild-type controls were set to 1 = 4). (D, F, H) Data are mean ± s.d. Source data for this figure is available on the online 
supplementary information page. 
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Figure 6 ChIP analysis of Sin3A/HDACl/HDAC2 target genes. Chromatin from littermate controls, Hdacl^^^^^ and Hdacl^^^^^ HdacI^^^^"^ 
epidermis was immunoprecipitated with antibodies specific for HDACl, HDAC2, Sin3A, histone H3K9ac, pan-histone H4ac or IgG as a control. 
Precipitated DNA was analysed by qRT-PCR with primers specific for Sprr2h, S100a9, Lce3b, Epgn and Vstm21 (control) promoter regions. Data 
are representative for at least two independent experiments. n^2-S. Mean value with s.d. is shown. 



compensate for the loss of HDACl in Hdacl^^^^^ mice, 
whereas in Hdacl^^^^"^ HdacZ^^^^"^ mice lower Sin3A and 
HDAC2 levels fail to fully exert the Sin3A/HDACl/HDAC2 
co-repressor function. 

Enhanced tumour development in K5-S0S Hdad^^^^^ 
mice 

HDACl and HDAC2 have been shown to be crucial for cell 
proliferation and are therefore considered as promising 
targets for anti-tumour therapy (Lagger et aU 2002; 
Haberland et aU 2010; Yamaguchi et aU 2010). However, 
the data shown above point towards an anti-proliferative 
function of HDACl during epidermal development and 
homeostasis. We therefore examined the effect of ablation 
of HDACl and HDAC2 in the genetic K5-S0S skin tumour 
model (Sibilia et aU 2000] . In this system, the ras pathway is 
constitutively activated in the epidermis due to a dominant 
form of Son of Sevenless (SOS) expressed under the control 
of the K5 promoter (K5-S0S] . In the presence of a functional 
epidermis growth factor receptor (EGFR), required to 
provide an essential survival signal to tumour cells, 
K5-S0S mice develop skin tumours that share features of 



human squamous cell carcinomas (Sibilia et al, 2000). 
Papillomas develop predominantly on the tail, the paw, 
behind the ears, and at sites subjected to scratching and 
biting. Notably, K5-S0S Hdacl^^^^^ mice showed accelerated 
onset of tumour development and a significant increase in 
relative tail tumour weight (Figures 7A-C). In accordance 
with our previous observations, deletion of HDACl led to an 
up-regulation of HDAC2 in K5-S0S Hdacl^^^^^ epidermis 
(Figure 7D). HDAC activity assays revealed a significant 
reduction in total cellular deacetylase activity in the epider- 
mis of the K5-S0S Hdacl^^^^"^ mice (Figure 7E). HDACl- 
deficient tumours showed hyperproliferation of the basal 
and suprabasal layers as determined by Ki67 staining and 
reduced expression of the epidermal differentiation marker 
KIO (Figures 7Fand G). Importantly, co-repressor associated 
deacetylase was significantly reduced for Sin3A, NuRD 
and Co-REST complexes in HDACl -deficient tumours com- 
pared to wild-type tumours (Figures 7H and I). In accor- 
dance with the link between Sin3A function and c-Myc 
expression, we detected increased levels of c-Myc protein 
and the c-Myc target Skp2 in KS-SOS Hdacl^^^^^ tumours 
(Figures 7J and K) . 
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Figure 7 Increased tumour development in K5-S0S Hdacl^^^^^ mice. (A) Picture of littermate control and K5-S0S Hdacl^^^^^ mice at Pll. 
Arrows indicate sites of premature tumour appearance on mutant mice. (B) Tumour appearance curves of control and K5-S0S Hdacl^^^^"^ mice 
(P = 0.0001) and (C) their relative tail tumour weight at Pll. ^ = 3-5. P = 0.0064. (D) Immunoblot analysis of HDACl and HDAC2 protein 
expression, with p-actin as a loading control and quantification of the immunoblot signal. 7i = 3. (E) HDAC activity measured in SOS tail 
epidermis of K5-S0S control and K5-S0S Hdacl^^^^^ mice. n = 4, P = 0.0128. (F) IHC staining and quantification of Ki67 in tail skin of control 
and K5-S0S Hdacl^^^^"^ (PH). Scale bar: SOjim. n^3, P = 0.0035. (G) Immunoblot analysis of KIO expression in epidermal extracts from 
control and K5-S0S Hdacl^^^^"^ mice, p-actin was used as a loading control. (H, I) Analysis of co-repressor complexes and associated 
deacetylase activity in control and K5-S0S Hdacl^^^^"^ mice. Input and immunoprecipitated proteins were analysed on immunoblots with 
HDACl, HDAC2, Sin3A, MTA2 and CoREST antibodies. p-Actin served as a loading control (H: ^ = 3, P = 0.0103, 0.003, 0.0358 and 0.0029). (J) 
Skin tumour sections of control and K5-S0S Hdacl^^^^^ mice (Pll) were labelled with c-Myc antibody and nuclei were counterstained with 
DAPI. Scale bar: 100 |im. (K) Expression of c-Myc and Skp2 in nuclear and whole-cell extracts respectively of control and K5-S0S Hdacl^^^^^ 
mice epidermal tumours. Lamin B and p-actin served as a loading control. An at least three-fold up-regulation of c-Myc was observed in three 
independent experiments. (C-F,H) Data represent means ± s.d. Source data for this figure is available on the online supplementary information 
page. 



In contrast, ablation of HDAC2 in the epidermis of 
K5-S0S mice had no effect on tumour appearance and tumour 
w^eight (Figures 8A-C). While HDACl protein levels were 



up-regulated in the absence of HDAC2, total cellular and 
co-repressor associated HDAC activities were unchanged in 
K5-S0S Hdac2^^^^^ epidermis (Figures 8D, E, H and I). 
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Figure 8 Loss of HDAC2 in the epidermis has no effect on tumour development in K5-S0S mice. (A) Picture of littermate control 
and K5-S0S HdacZ^^^^^ mice at Pll. (B) Tumour appearance curves of control and K5-S0S Hdac2^^^^^ mice. Note that in the genetic 
background of Hdac2^^^ mice tumours appear earlier compared to background of Hdacl^^^ mice. P= 0.8829. (C) Relative tail tumour weight of 
control and K5-S0S Hdac2^^^^^ mice at Pll. n = 3-5, P = 0.7293. (D) Immunoblot analysis of HDACl and HDAC2 protein expression, with 
P-actin as a loading control and the quantified immunoblot signal. n = 3. (E) HDAC activity measured in SOS tail epidermis of K5-S0S control 
and K5-S0S HdacZ^^^^^ mice, n = 3, P = 0.3996. (F) IHC staining and quantification of Ki67 in tail skin of control and K5-S0S Hdac2^^^^^ (Pll). 
Scale bar: 50|im. n = 3, P = 0.8483. (G) Immunoblot analysis of KIO expression in protein extracts from control and K5-S0S Hdac2^^^^^ 
epidermis. |3-actin was used as a loading control. (H, I) Analysis of co-repressor complexes and associated deacetylase activity in control 
and K5-S0S Hdac2^^^^'^ mice. Input and immunoprecipitated proteins were analysed on immunoblots with HDACl, HDAC2, Sin3A, MTA2 
and CoREST antibodies. |3-actin served as a loading control. (H: n = 3, P = 0.2697, 0.6416, 0.9631 and 0.4477). (J) Skin tumour sections of 
control and K5-S0S Hdac2^^^^'^ mice (Pll) were labelled with c-Myc antibody and nuclei were counterstained with DAPI. Scale bar: 100 |im. 
(K) Expression of c-Myc and Skp2 in nuclear or whole-cell epidermal extracts respectively of control and K5-S0S Hdac2^^^^^ mice tumours. 
Lamin B and p-actin served as a loading control. (C-F, H) Mean value with s.d. is depicted. Source data for this figure is available on the online 
supplementary information page. 



In agreement with the unaffected tumour size, expression of 
c-Myc, Skp2 and the proliferation marker Ki67 was 
comparable to littermate controls whereas the differentiation 
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marker KIO was slightly reduced in K5-S0S Hdac2^^^^^ 
tumours (Figures 8F, G, J and K). Thus, HDACl but not 
HDAC2 negatively controls tumour cell proliferation in the 
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epidermis. These results point towards a novel and unexpected 
role of HDACl as a tumour suppressor in the epidermis. 

Discussion 

Individual and redundant functions of HDACl and 
HDAC2 in the epidermis 

In this report, we analysed individual and redundant func- 
tions of HDACl and HDAC2 in epidermal development, 
homeostasis and tumorigenesis. Strikingly, K5-Cre mediated 
deletion of three of the four Hdacl /Hdac2 alleles had different 
consequences depending on the remaining allele. Mice with a 
single Hdacl allele in the epidermis [Hdacl^^^^^^ Hdac2^^^^^) 
showed no significant changes in HDAC activity, proliferation 
and differentiation in the epidermis, indicating that one 
Hdacl allele is sufficient to maintain proper epidermal 
development. In contrast, mice with a single HdacZ allele 
{Hdacl^^^^"^ HdacZ^^^^^] displayed a severe developmental 
phenotype accompanied by a significant reduction in total 
HDAC activity. These data reveal a predominant role of 
HDACl during epidermal development. The severe develop- 
mental phenotype of Hdacl^^^^^ Hdac2^^^^^ mice using 
K5-Cre was also described in a recent report using K14-Cre 
albeit without description of a molecular mechanism 
(Hughes et al, 2013] but were not uncovered in another 
report using the same Cre transgene (LeBoeuf et al, 2011). 
The reasons for this discrepancy are unclear. 

HDACl was identified as the major deacetylase during 
embryonic development (Lagger et al, 2002], in B cells 
(Reichert et al, 2012) and T cells (Grausenburger et al, 2010; 
Dovey et al, 2013; Heideman et al, 2013). On the contrary, a 
single Hdac2 but not Hdacl allele has been shown to be 
sufficient for normal oocyte (Ma et al, 2012) and brain 
development (Hagelkruys, Lagger et al, manuscript in 
revision). Taken together, these findings indicate that 
HDACl and HDAC2 do not only have redundant but also 
specific functions in differentiation and development. 

Loss of HDACl in the absence of compensating increased 
levels of HDAC2 in Hdacl^^^^^ Hdac2^^"^^^ epidermis led to a 
severe reduction in HDAC activity. This was in part due to 
destabilization of certain components of the Sin3 and NuRD 
co-repressor complexes such as Sin3A, MTAl and MTA2. 
Thus, in addition to their catalytic functions HDACl and 
HDAC2 seem to also have a scaffolding function required 
for the stability of epidermal co-repressor complexes. 
Notably, a similar dosage-dependent effect of HDACl/ 
HDAC2 ablation on HDAC activity and repressor complex 
function was recently observed in mouse T cells (Dovey et al, 
2013). HDACl -deficient T cells with only one Hdac2 allele 
exhibited a significant reduction in HDAC activity associated 
with Sin3A and NuRD complexes and reduced Sin3A and 
MTA2 protein levels. These findings support the hypothesis 
that beyond their enzymatic function HDACl/2 are also 
crucial for the structural integrity of repressor complexes. 

Control of epidermal proliferation, differentiation and 
lineage faith by HDACl/2 co-repressor complexes 

Hdacl^^^^^Hdac2^^^^'^ mice exhibited profound changes 
in the development of all three epidermal lineages resulting 
in epidermal hyperproliferation and hyperkeratosis, accom- 
panied by hair loss and SO enlargement. The most remark- 
able phenotype of mutant mice is the dramatic alopecia. 
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Hdaci^/^^P Hdac2^^^^'^ HFs are generated with similar 
frequency but fail to properly enter and transit the first hair 
cycle. Already during morphogenesis mutant HFs show 
increased levels of p53 and apoptosis resulting in HP atrophy. 
Accordingly, the absence of HDAC1/HDAC2 in the epidermis 
was shown to lead to high levels of acetylated p53 (LeBoeuf 
et al, 2011) suggesting a direct or indirect role of the two 
class I deacetylases as negative regulators of p53 activity. 
Interestingly, in contrast to mutant HFs the IFE of Hdacl^^^^"^ 
Hdac2^^^^^ mice is more resistant to p53 induced apoptosis. 
This might be one explanation for the hneage-specific 
differences in the Hdacl^^ ^^^Hdac2^^ ^ epidermis. 

On the other hand, up-regulation of c-Myc, a central 
regulator of epidermal lineage commitment, contributes to 
impaired epidermal homeostasis of Hdacl^^^^^ Hdac2^^^^^ 
mice. Overexpression of c-Myc was shown to mobilize SCs 
from the SO compartment and to promote proliferation as 
well as differentiation of the IFE and SO lineage at the 
expense of HF development (Arnold and Watt, 2001; 
Lo Celso et al, 2008; Berta et al, 2010). Accordingly, we 
observed SO hypertrophy and increased levels of the SO 
marker Lgr6 in mutant mice. We cannot fully exclude that 
wounding contributes to the mobilization of HF SCs to the 
IFE of Hdacl^^^^^ Hdac2^^^^'^ mice; however, we could not 
detect increased numbers of infiltrating immune cells in the 
mutant epidermis. Hdacl^^^^^Hdac2^^ mice phenocopy to 
a large extent mice with epidermal deletion of Sin3A 
(Nascimento et al, 2011). We propose that the partially 
impaired function of the Sin3A co-repressor complex in the 
Hdacl^^^^^Hdac2^^^^^ epidermis has two consequences. 

First, this leads to stabilization and increased expression of 
the c-Myc protein. Sin3A was shown to counteract c-Myc 
activity by causing the deacetylation of c-Myc protein 
(Nascimento et al, 2011), while acetylation has been shown 
to enhance the stabihty of the c-Myc protein (Vervoorts et al, 
2003; Patel et al, 2004). In agreement with the idea that the 
Sin3A/HDACl/HDAC2 co-repressor complex attenuates 
c-Myc activity by reducing protein stabihty, we found that 
the c-Myc protein can be stabilized by HDAC inhibition and 
was increased in the Hdacl^^^^^Hdac2^^ epidermis 
without changes in mRNA levels. In Hne with an important 
role for c-Myc as a contributor to the impaired epidermal 
development, we discovered a significant overlap between 
genes deregulated in Hdacl^^^^^Hdac2^^ epidermis with 
genes deregulated in c-Myc overexpressing mice (Frye et al, 
2003). 

Second, the reduced amount of Sin3A and the impaired 
Sin3A co-repressor function led to derepression of differentia- 
tion specific genes, thereby contributing to the hyperkeratosis 
of the mutant epidermis. ChIP experiments show that 
reduced Sin3A recruitment and uncompensated loss of 
HDACl at EDO genes resulted in increased local histone 
acetylation and gene expression. Given the partially impaired 
function of CoREST and NuRD complexes and the large 
number of deregulated genes other HDAC1/HDAC2 recruiting 
factors and complexes in addition to the Sin3A complex 
might also play important roles in epidermal development 
and homeostasis. 

HDACl acts as a tumour suppressor in the epidermis 

In an SOS-driven skin tumour model deletion of HDACl in 
the epidermis led to enhanced tumorigenesis associated with 
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increased proliferation and reduced differentiation, whereas 
ablation of HDAC2 had no effect. Hdacl^^^^^ mice did 
not display an obvious developmental phenotype but the 
frequently observed scar formation at the tail and the accel- 
erated tumour development in the K5-S0S tumour model 
suggest that under mechanical or oncogenic stress conditions 
HDAC2 cannot fully compensate for the loss of HDACl in the 
epidermis. In accordance with this idea, total deacetylase and 
co-repressor activity was unchanged in Hdacl^^^^^ epidermis 
but clearly reduced in HDACl -deficient tumours. In particu- 
lar, the impaired activity of the Sin3A co-repressor complex 
and the concomitant overexpression of the c-Myc in HDACl- 
but not HDAC2- deficient papillomas provide strong evidence 
for crucial function for Sin3A and HDACl as a negative 
regulator of the proto-oncogene c-Myc. Thus, the hyperpro- 
liferation and spontaneous tumour development observed in 
Hdacl^^^^'^HdacZ^^^^'^ epidermis and in Hdaci -deficient SOS 
skin tumours indicates a novel role of HDACl as a tumour 
suppressor in the skin. These observations show parallels 
with two recent studies on the function of HDAC1/HDAC2 in 
T cells. Deletion of two Hdacl alleles and a single Hdac2 allele 
by Lck-Cre resulted in neoplastic transformation of immature 
T cells (Dovey et al, 2013; Heideman et aU 2013]. In this study, 
enhanced c-Myc expression was detectable at the mRNA level 
and mostly due to trisomy of chromosome 15. The authors 
concluded that HDACl and HDAC2 control genomic stabihty 
(Dovey et al, 2013; Heideman et al, 2013] and p53 activity 
controlling functions (Heideman et al, 2013] in a dosage- 
dependent manner. 

The function of HDACl as a negative regulator of prolifera- 
tion in the epidermis is in contrast to some earlier studies in 
which HDACl deficiency in ES cells or fibroblasts showed 
reduced rates of proliferation due to increased levels of the 
CDK inhibitor p21 (Lagger et al, 2002; Zupkovitz et al, 2010]. 
Similarly, loss of HDACl in human tumour cells led to 
reduced proliferation and elevated apoptosis (Glaser et al, 
2003; Senese et al, 2007]. On the contrary, it was shown that 
the conditional loss of HDACl in T cells results in increased 
proliferation (Grausenburger et al, 2010] and HDACl - 
deficient murine teratomas display enhanced proliferation 
and reduced differentiation (Lagger et al, 2010] . 

In summary, we have discovered a crucial function of 
HDACl in epidermal homeostasis and as a tumour suppressor 
in skin cancer. Our results have potential consequences for 
medical applications and the use of HDAC inhibitors in 
tumour therapy. For instance, low expression of HDACl/ 
HDAC2 in the epidermis might cause a predisposition to- 
wards neoplasia. Furthermore, HDAC inhibitors that mostly 
target HDACl could induce undesired effects such as hyper- 
proHferation. Finally, drugs that target both HDACl and 
HDAC2 might be promising since loss of both enzymes 
invariably induces cell death in all proliferating cell types 
and tissues tested so far. 

Materials and methods 

Animal care and transgenic mouse lines 

All mouse lines were bred to a mixed genetic background of C57BL/ 
6J X 129SV. Keratin 5 Ore [K5-Cre) mice (Ramirez et al, 2004] were 
mated to Hdacl^^^ and Hdac2^^^ mice (Yamaguchi et al, 2010] to 
generate mice with deletions of Hdacl and/or Hdac2 alleles in the 
epidermis. K5-S0S tumour mice with hypomorphic EGFR alleles 
^£(.pj^wa2/wa2^ (Sibilia et al, 2000] were crossed to Hdacl^^^ and 
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Hdac2^^^ mice. Mice were kept in the animal facilities of the MFPL 
and the Medical University of Vienna in accordance with 
institutional policies and federal guidelines. 

Epidermis isolation and keratinocyte culture 

Murine tail skin was isolated and incubated dermis side down in 
PBS containing 5 mg/ml Dispase (Roche] for 2 h at 4°C. Epidermis 
was removed from dermis, washed in PBS and snap frozen. Tissue 
was further used for protein or RNA isolation. Mouse epidermal 
keratinocytes were isolated as previously described (Sibilia et al, 
2000] and cultured in keratinocyte growth medium (KGM 
BulletKit®, Clonetics®) containing 8% chelated fetal calf serum 
and 0.05 mM CaCl2 on collagen-fibronectin coated dishes. 

RNA isolation and real-time PCR analysis 

Total RNA was isolated following the manufacturer's instructions 
(TRIzol, Invitrogen]. RNA was reversely transcribed with the iScript 
cDNA synthesis kit (Bio-Rad]. Real-time PCR analysis was per- 
formed with the KAPA SYBR FAST qPCR MasterMix (Peqlab] on 
an iCycler IQ system (Bio-Rad]. Relative expression levels were 
normalized to HPRT. For primers, see Supplementary data. 

Agilent microarray data and overlaps 

Gene expression profiling was performed with Agilent Whole 
Genome Microarrays (Agilent's SurePrint G3 Mouse GE 8 x 60K 
Microarray, Agilent Microarray Design ID 028005]. The data analy- 
sis was performed using GeneSpring software 11.5 (Agilent 
Technologies]. The cutoffs for differential expression were set as 
absolute fold-change >2 and a corrected P-value of <0.05. GO 
analysis was performed with the DAVID (http://david.abcc.ncifcrf.- 
gov/]. Overlaps were identified with R and compared to the mean of 
1000 random overlaps. The significance (P-value] from these results 
was calculated using Fisher's exact test, a hypergeometric prob- 
ability test to determine the significance between the overlaps 
(R version 2.15.0 (2012-03-30]; function fisher.test]. 

Protein isolation, immunoblot analysis and HDAC activity 
assays 

For protein isolation from cultured keratinocytes, cells were scraped 
and isolated in Frackelton buffer (10 mM Tris-HCl pH7.05, 50 mM 
NaCl, 30 mM sodium pyrophosphate, 50 mM NaF, 1% Triton X-100]. 
Protein extracts were vortexed and incubated on ice for 10 min. For 
protein isolation from tissue, the epidermis was isolated as described 
above. The tissues were homogenized in Frackelton buffer with 
ceramic beads (Precellys] for 30 s at 6m/s with a FastPrep-24 
homogenizer (MP Biomedicals]. Equal amounts of protein extracts 
were separated by SDS-PAGE and transferred onto nitrocellulose 
membranes (Protran, Whatman] according to the standard protocols. 
For detection, the enhanced chemiluminescence kit (Perkin-Elmer] 
was used. HDAC activity assays were performed with epidermal 
protein extracts as previously described (Lagger et al, 2002]. 
Frackelton buffer was supplemented with protease inhibitor 
cocktail (complete, Roche] and aprotinin (Sigma]. For antibodies, 
see Supplementary data. 

Chip and PCR analysis 

Epidermis was chopped, washed with PBS and crosslinked with 
disuccinimidyl glutarate (2mM, AppliChem] for 25 min at room 
temperature (RT]. After a PBS washing step, the tissue was cross- 
linked with formaldehyde (final cone, of 1%] at RT for 10 min. 
Crosslinking was stopped by addition of 125 mM glycine. Chromatin 
isolation was done as previously described (Hauser et al, 2002]. For 
Chip, equal amounts of sonicated chromatin were diluted 10-fold 
and precipitated overnight with antibodies listed in Supplementary 
data. Chromatin antibody complexes were isolated using protein 
A-beads or G-beads (Dynabeads®, Invitrogen] and the extracted 
DNA was used for qRT-PCR. In parallel, PCRs with 1:20 dilutions 
of genomic DNA (input] were carried out. For primers, see 
Supplementary data. 

Co-immunoprecipitation assay 

Total protein extracts from keratinocytes or epidermal extracts were 
harvested in Frackelton lysis buffer. Equal amounts of protein were 
incubated for 2 h at 4°C with 4 |ig antibody. Immunoprecipitation 
was carried out using protein A-beads or protein G-beads 0/N at 
4°C. Immune complexes were washed with Frackelton lysis buffer. 
Samples were used for HDAC activity assays or immunoblotting. 
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Nuclear extraction 

Isolated epidermis was chopped and filtered through a 70 \xm cell 
strainer. After washing, the cell pellet was resuspended in sucrose 
buffer (0.32 M sucrose, 10 mM Tris-HCl pH 8.0 3mM CaCl2, 2mM 
MgOAc, 0.1 mM EDTA, 0.5% NP-40). Nuclei were isolated by 
centrifugation and the nuclear pellet was washed twice with sucrose 
buffer (without NP-40). Then, the pellet was resuspended in low 
salt buffer (20 mM HEPES pH 7.9, 1.5 mM MgCl2, 20 mM KCl, 
0.2 mM EDTA, 25% glycerol) and mixed with an equal amount of 
high salt buffer (1.5 mM MgCl2, 800 mM KCl, 0.2 mM EDTA, 25% 
glycerol, 1% NP-40). After shaking for 1.5 h at 4°C the nuclear 
fraction was obtained as a supernatant after centrifugation. 

Histological and IHC analyses 

Tissue samples were fixed overnight in 4% PFA and further 
embedded in paraffin. Stainings were performed on 4 |im sections. 
H&E stainings were carried out according to the standard procedure 
with an ASSl staining unit (Pathisto). IF stainings were carried out 
with the Tyramide Signal Amplification Kit (PerkinElmer) according 
to the manufacturer's instruction. The slides were mounted 
with DAPl in Vectashield (Vector Laboratories). For IHC stainings 
with mouse primary antibodies the M.O.M Immunodetection Kit 
(Vector Laboratories) and for rabbit primary antibodies the 
VECTASTAIN ABC Kit (Vector Laboratories) were used according 
to the manufacturer's protocol. Antibody detection was carried out 
with the DAB Peroxidase Substrate Kit (Vector Laboratories) and 
sections were counterstained with haematoxylin reagent. For Oil O 
red staining of SG and for quantification of histology stainings, see 
Supplementary data. 

Microscopy 

IHC stainings, tissue whole mounts and hair samples were captured 
on a Zeiss stereomicroscope with camera. IHC fluorescence stain- 
ings were imaged on a Zeiss LSM Meta 510 confocal microscope. 

Statistical analysis 

qRT-PCR and ChIP experiments were evaluated with Microsoft 
Excel. Immunoblot signal intensities were quantified using the 
ImageQuant® software and relative protein levels were normalized 
to p-actin or Lamin B. The significance between groups was 
determined by the unpaired Student's t-test. P-values were calcu- 
lated with the Graph-Pad Prism software and standard deviation 
(s.d.) is shown. *P<0.05; **P<0.01; ***P<0.001; ns = not signifi- 
cant. In Figures 7B and 8B, the data were analysed by a log-rank 



(Mantel-cox) test. If data were obtained from independent experi- 
ments, then controls were set to 1 or 100% and the s.d. of the ratio 
of the mutant mice in relation to the controls is shown. 

For further information of materials (primers, antibodies) and 
methods, see Supplementary data. 

Supplementary data 

Supplementary data are available at The EMBO Journal Online 
(http://www.embojournal.org) . 
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